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Introduction 


\ 

In  the  last  five  years,  much  progress  has  been  made  in  the  study  of 

acoustic  wave  scattering  by  elastic*  obstacles  immersed  in  water  using  the 

T-matrix  approach.  Excellent  agreement  with  experiments  was  demonstrated 

\  12 

for  the  finite  elastic  cylinder  ’  with  spherical  end  caps.  The  T-matrix 

3 

approach  was  then  extended  to  the  problem  of  elastic  shells  immersed  in 

V, 

water.  Numerical  calculations  of  the  frequency  dependence  of  acoustic 
wave  scattering  by  prolate  and  oblate  sphero^da'l  elastic  and  viscoelastic 
shells  were  presented  in  Ref.  3.  However,  since  the  problem  was  formulated 
exactly  in  the  shell  region,  several  numerical  difficulties  arose  as  the 
thickness  of  the  shell  was  decreased  and  the  frequency  increased.  .It  is 
understood,  that  in  problems  of  interest  to  the  Navy,  scattered  field  data 
is  required  for  long,  thin  bodies  of  revolution  whose  wall  thickness  is 
very  small  compared  to  the  wavelength  of  the  incident  wave  but  whose  overall 
dimensions  are  comparable  and  often  larger  than  the  wavelength  of  the  incident 
wave. 

This  project  was  begun  to  precisely  address  such  problems,  taking  into 
account  and  exploiting  the  thinness  of  the  shell.  We  propose  to  use  shell 
theory  equations  rather  than  the  full  elasticity  equations  in  region  I  (see 
figure  1)  in  order  to  avoid  some  of  the  complications* encountered  in  Ref.  3. 
Shell  theory  reduces  the  partir  *^erential  equations  in  three  coordinates 
to  a  higher  order  equation  in  two  coordinates.  All  description  of  the 
displacement  and  stress  fields  is  made  with  respect  to  a  reference  surface 
in  the  shell.  Now,  it  is  no  longer  necessary  to  invoke  integral  repre¬ 
sentations  to  describe  the  elastic  field  in  this  region.  One  extracts 
from  the  shell  equations  an  expression  for  the  impedance  matrix  of  the 
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shell  which  is  then  interfaced  with  the  T-matrix  approach  for  acoustic 
wave  scattering  problems.^: - 


The  null  field  equations 

Consider  a  closed  obstacle  of  piecewise  smooth  surface  S  with  an  out- 

A 

ward  normal  n  which  is  immersed  in  water.  A  plane  harmonic  acoustic  wave 

of  frequency  u  is  incident  on  the  obstacle.  The  total  field  in  the  water 

consists  of  the  incident  field  <p^  and  the  scattered  field  The  integral 
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representation  of  these  fields  is  well  known  ’  and  takes  the  form 


-(p  (r);  r  inside  S 


i^gCr);  r  outside  S 


V'g(r,r’))()»_^(r') 


g(r.r')  (n'  ♦  V 


dS(r’) 


where  g(r,r')  is  the  Green’s  function  of  the  scalar  Helmholtz  equation 
for  an  infinite  medium  and  is  the  velocity  potential.  The  expressions 
for  'g'  are  known  in  both  2-D  and  3-D 

g(r,r')  a  iiT  (klr-r’|)  ;  2-D 


g(r,r')  a  exp  i  k)r-r'|  /  kjr-r'  ;  3-D 


In  the  T-matrix  approach  we  expand  the  incident  and  scattered 
fields  in  eigen  functions  of  the  wave  equations  as  follows 


(|»^(r)  a  I  A^Re<>^(r) 
n 


♦j(r)  a  I  fjj0u4>jj(r) 


*  ro;'Y 

TtOj 


t,: 

«3  «  I 

5  g  f  s  i  J 

5  ^ fj  w  j 


«i  o 
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where  <J(^  are  the  circular  cylindrical  functions  in  2-D  and  the  spherical 
functions  in  3-D.  The  index  'n'  is  a  super  index  that  represents  all 
the  required  subscripts  for  the  wayefunction.  The  qualifier  ’Re'  stands 
for  functions  that  are  regular  at  the  origin  and  'Ou*  represents  functions 
that  are  outgoing  at  infinity.  Thus  these  expansions  are  consistent  with 
our  knowledge  of  the  regularity  of  the  incident  field  in  the  region  occupied 
by  the  obstacle  and  the  radiation  conditions  satisfied  by  the  scattered 
field  at  distances  far  from  s.  A  detailed  description  of  the  wavefunctions 
for  2-D  and  3-D  problems  may  be  found  in  Ref.  4.  It  may  be  noted  that 
since  only  scalar  waves  are  involved,  we  may  use  elliptical  or  spheroidal 
functions  as  appropriate  for  the  e:q>ansion  of  the  fields^.  Thus  repre¬ 
sents  any  complete,  orthogonal  set  of  solutions  of  the  scalar  Helmholtz 
equation. 

The  expansion  of  the  Green's  function  is  also  known  and  takes  the  form 


gCr,r»)  =  i  ^  (Xi  Re 


(6) 


where  r^  and  r^  stand  for  the  greater  and  lesser  of  r  and  r'  respectively. 
The  factor  tt  is  used  in  2-D  and  the  factor  k  in  3-D 

Substituting  eqs.  (4) -(6)  in  eq.  (1)  we  arrive  in  the  usual  manner  at 


(7) 


In  order  to  solve  for  the  known  scattered  field  coefficients  'f^' 
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in  terms  of  the  knovm  incident  field  coefficients  we  must  specify 
and  n  •  on  the  surface  S  of  the  elastic  shell.  Since  41  is  the 
velocity  potential,  for  time  harmonic  waves,  is  related  to  the  surface 
pressure,  and  n  •  7^41  to  the  normal  component  of  the  particle  velocity, 
n  •  v^  on  the  surface.  Thus,  since 

p(r)  =  i  u)  p  4>  (r)  (8) 


and 


v(r)  =  V4»(r) 


C9) 


equation  (7)  may  be  rewritten  as 


IB*  ■  • 

i ■  =-0 


dS 


where  p  is  the  mass  density  of  the  water.  We  now  proceed  to  discuss  the 
use  of  shell  theory  to  represent  the  pressure  and  normal  velocity  on  the 
surface  of  the  shell  keeping  in  mind  that  the  boundary  conditions  at  S 
are 


(10) 


P^(r)  a  P^  (r)  ;  r  on  S 


(11) 


and 


n  •  v^(r)  a  n  .  v_(r)  ;  r  on  S. 


(12) 


Impedance  matrix  for  elastic  shells 


Consider  a  coordinate  system  on  the  surface  S  that  is  2-D  and  tangent 
to  the  surface  at  every  point.  Let  C5>ri)  represent  these  two  coordinates. 

Let  u^  be  the  components  of  the  displacement  field  in  C,  5  and  n  directions 
respectively,  where  c  is  a  coordinate  normal  to  the  shell  surface.  If  the 
thickness  '2h'  of  the  shell  is  small  compared  to  its  other  dimensions,  it 
is  appropriate  to  assume  that  u,  v  and  w  are  functions  of  S  and  n  only.  Let 
i(i^  be  a  complete  set  of  functions  that  are  orthogonal  on  the  surface  S. 

Then  we  may  write 

“its-  3  *  I  )  ;  i  =  1,2,3  (13) 


The  set  of  coefficients  can  be  considered  the  generalized  coordinates 
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of  the  problem.  This  approach  has  been  previously  used  by  Junger  to 
analyze  the  scattering  of  acoustic  wave  circular  cylindrical  shells  of 
infinite  length  when  waves  are  incident  normal  to  the  axis  of  the  cylinder. 
We  may  set  up  the  Lagrangian  for  the  shell  as 

L  =  T  -  V 

where  T  is  the  kinetic  energy  density  and  V  the  potential  energy  density, 
due  to  the  elastic  strain  energy  of  the  shell.  Hence 

T  *  hp  I  (Uj  ♦  U2  ♦  UjJdS 
S 


(14) 


(15) 


(16) 
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The  function  F  that  is  quadratic  in  the  displacements  may  be 
constructed  from  shell  theory. 

In  shell  theory  approximations,  one  considers  various  contributions 
to  the  elastic  strain  energy  that  includes  effects  such  as  the  membrane 
effect,  bending,  rotary  inertia,  transverse  shear  etc.  These  effects  must 
be  included  as  appropriate  for  the  particular  problem  on  hand  with  con¬ 
sideration  given  to  the  frequency  of  the  incident  harmonic  wave  and  the 
exact  shape  of  the  shell.  At  the  present  time  this  part  of  the  study  is 
still  incomplete.  However,  this  does  not  prevent  us  from  proceeding  further 
with  the  formulation. 

Lagrange's  equations  for  the  system  takes  the  form 


(17) 


(18) 


(19) 


Where  is  the  generalized  force  due  to  the  pressure  exerted  by  the 
fluid  outside  the  shell.  Since  the  fluid  is  assumed  to  be  non-viscous, 
the  force  acts  only  in  the  normal  direction  i.e.,  along  the  coordinate  c* 

The  generalized  force  in  the  other  two  directions  are  zero. 

An  expression  for  can  be  obtained  by  giving,  the  shell  a  vertical  dis¬ 
placement  8u  in  the  ;  direction.  If  is  the  total  pressure  on  the  shell 
surface,  then 


Qn  = 


(20) 


The  total  pressure  on  S  is  expanded  as  follows 


=  I  J  ^OnS 

m 


(21) 


Thus 


On  “  ■  ^  j  '('mC^*n)«^„(5.n)dS 

m  i 


(22) 


If  the  4i  are  an  orthonozmal  set 
n 


On 


(23) 


We  further  note  using  eq.  (13)  that 


(24) 


and 


3V 

3o^ 


3  f  ..  3  .. 

»  y  u.t|»^dS  *  I  a_ 

j^l  i  n  1  n  n  n 


(25) 


where  the  3x3  matrix  A^^  may  be  set  up  for  each  modal  index  n  from  the 
explicit  expressions  for  the  strain  energy  density  of  the  shell.  We 
repeat  that  this  part  of  the  formulation  is  inconplete. 

Substituting  eqs.  (24)  and  (25)  in  Lagrange's  equations  as  given  in 
eqs.  (17)  -  (19),  we  obtain 
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Ij  j  - 


•2hPsW  “n  +  I 


2j  J  = 


.2hPga)  +  I  A^  CL^ 


-2h^/4  "  I  “i  “  ° 

S  n  ^ 


Equations  (26)  -  (28)  may 


be  solved  to  yield  the  ratio  which  we 


write  as 


3:  .» 

a  1  n 
n 


from  eqs.  (13)  and  (21)  we  have 


Further  from  the 


boundary  conditions  on  S  as  given  in  eq.  (12) ,  we  have 


n  •  %  =  It  "l^^^  ’  ^  ^ 

Again  using  eq.  (13) 

i  .  V  =  -  iu  I 
n 

The  medimical  of  the  shell  cm  b.  defined  as 


n  •  v. 
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and  from  eqs.  (30)  and  (32)  -  (33)  we  can  define  a  modal  impedance  Z^. 
Equations  (29)  and  (34)  may  now  be  used  in  the  null  field  equations  (10) 
to  relate  the  scattered  and  incident  field  coefficients. 

The  T -matrix  of  a  thin  elastic  shell 

Using  eqs.  (30)  and  (32)  in  (10),  we  have 


where 


Z  =  -iz  /doj 
m  n 


Equation  (34)  may  now  be  solved  in  the  usual  manner  to  yield 


n 


T  a 
nm  m 


(35) 


(36) 


where 


T  =  Q(Q)'^ 
d 

1^1’ 
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Closure 


At  this  time  the  formulation  of  the  scattering  problem  using  impedance 
boundary  conditions  is  complete.  However,  the  exact  form  of  the  nodal 
impedance  must  still  be  extracted  from  the  shell  equations  for  the  strain 
energy  density.  After  this  has  been  done  for  shells  of  various  shapes, 
numerical  calculations  must  be  performed  for  various  model  shapes,  so  that 
they  can  be  compared  with  available  exact  calculations  (using  the  full 
elasticity  equations). 
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